Gap solitons under competing local and nonlocal nonlinearities 
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We analyze the existence, bifurcations, and shape transformations of one-dimensional gap solitons 
(GSs) in the first finite bandgap induced by a periodic potential built into materials with local self- 
focusing and nonlocal self-defocusing nonlinearities. Originally stable on-site GS modes become 
unstable near the upper edge of the bandgap with the introduction of the nonlocal self-defocusing 
nonlinearity with a small nonlocality radius . Unstable off-site GSs bifurcate into a new branch 
featuring single-humped, double-humped, and fiat-top modes due to the competition between local 
and nonlocal nonlinearities. The mechanism underlying the complex bifurcation pattern and cutoff 
effects (termination of some bifurcation branches) is illustrated in terms of the shape transformation 
under the action of the varying degree of the nonlocality. The results of this work suggest a possibility 
of optical-signal processing by means of the competing nonlocal and local nonlinearities. 
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I. INTRODUCTION 

The concept of photonic crystals (PhCs), i.e., artificial 
media with a periodic modulation of local optical charac- 
teristics, provides ways to tailor the dispersion, diffrac- 
tion, and routing of electromagnetic waves [1^. As for 
natural crystals, fundamental characteristics of the PhCs 
are the band-diagrams, which reveal gaps where Bloch 
waves cannot propagate. In PhCs made of nonlinear ma- 
terials, self-trapped localized modes in the form of the 
gap solitons (GSs) may exist, as a result of the interpla" 
of the Kerr-type nonlinearity and periodic structures 

Unlike spatial bright solitons supported by the bal- 
ance between the self-focusing nonlinearity and diffrac- 
tion in uniform bulk media [6], the dispersion relation 
induced by the PhC makes it possible to create GSs in 
both focusing and defocusing media. Combining assets 
of PhCs and regular solitons, GSs have a potential for ap- 
plications to soliton-driven photonics. New technologies 
enabling reconfigurable optical lattices, such as photore- 
fractive crystals f?] and nematic liquid crystals [sj, also 
open new ways to control the dynamics of solitary waves 
by adjusting the lattice depth and period. 

While the modulational ^ and oscillatory instabilities 
[Tol [TTI | impose limits on the use of the GSs, the stability 
and mobility of the GSs may be enhanced in nonlinear 
media featuring a nonlocal response (l^ - [T5| . The nonlo- 
cal nonlinearity is important when the correlation radius 
of the material's response function becomes comparable 
to the transverse width of the wave packet [l^ ■ The 
nonlocal nonlinearity gives rise to specific features in the 
soliton dynamics, including the modification of the mod- 
ulational [181], azimuthal ^QJ, and transverse |20|| insta- 
bilities. Sup pression of the collapse of multidimensional 
solitons |2l| , a change of the character of interactions be- 
tween them [2^, formation of soliton bound states [23j . 
merger of colliding solitons into a standing wave [13] , and 
families of dark-bright soliton pairs (25j were predicted 



recently too. Experimental observations of the nonlocal 
response have also been demonstrated in sundry media, 
including photorefractive media [2^ , nematic liquid crys- 
tals [l?!, and thermo-optical materials [1^, [l^j, with a 
large tunable range of the nonlocality degree. 

In this work, we aim to study GS modes in the first 
bandgap of the model including local self-focusing and 
nonlocal self-defocusing nonlinearities. When the non- 
locality radius is zero, we assume equal magnitudes of 
the self- focusing and self-defocusing terms, i.e., complete 
cancelation of the nonlinearity, hence no solutions in the 
bandgap. For a small degree of the nonlocality, the exis- 
tence, bifurcation, and shape transitions of the emerging 
bright GSs are analyzed. With the competing local and 
nonlocal interactions of opposite signs, the family of on- 
site GS modes remain stable, obeying the " anti-Vakhitov- 
Kolokolov" criterion js^], for the case that the nonlocal 
perturbation of the refractive index is small, but become 
unstable near the upper edge of the bandgap. However, 
the bifurcation generating the GSs near the other edge of 
the bandgap features an inverted slope of the bifurcation 
curve for a relatively small degree of the nonlocality, in 
which case the on-site GSs are unstable. The off-site GS 
family features single-humped, double-humped, and flat- 
top profiles for different degrees of the nonlocality. We 
also investigate the situation in the space of the soliton's 
propagation constant and power, for the varying nonlo- 
cality degree, in order to illustrate the mechanism under- 
lying the complex bifurcation pattern and related cutoff 
effects (termination of some solution branches). Using 
results reported in this work, we discuss the possibility 
to design GS-based signal-processing schemes by dint of 
manipulating the nonlocal interactions. 

In addition to optics, GSs of matter waves have also 
been theoretically studied 
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and experimentally cre- 
ated [3^ in Bose-Einstein condensates (BECs) formed 
by atoms with repulsive interactions, trapped in optical- 
lattice (OL) potentials. In addition to the known contact 
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interaction in the BECs of alkali-metal atoms, the inter- 
action of chromium atoms, ^^Cr, includes a dipole-dipole 
interaction, which is intrinsically anisotropic and nonlo- 
cal. The condensate of ^^Cr was created and investigated 
using magnetic (33i-i36i] and all-optical [STjj traps, see also 
review [38|. By adjusting the orientation of the dipoles, 
one can effectively control the nonlocal dipole-dipole in- 
teractions. For the dipolar BEC trapped in OLs, the 
competition between the contact and long-ranged dipole- 
dipole interactions not only dramatically change the band 
structures of nonlinear Bloch waves [391 , but also modi- 
fies families of matter- wave solitons |40|. 

The rest of the paper is organized as follows. In SecHIl 
the model including the competing local and nonlocal 
nonlinearitics is described. Properties of the on-site and 
ofF-site GS families supported by the local self-focusing 
nonlinearity are recapitulated to show a transition of 
mode profiles in the first bandgap. Results produced 
by the interplay of the local self-focusing and nonlocal 
self-defocusing nonlinearitics for on-site and ofF-site GS 
families are reported in Sections IIIII and ITVl respectively. 
Tracing the change of the corresponding GS shapes in 
the parameter planes, we explain the character of the 
corresponding bifurcation and identify a possible control 
mechanism for the optical-signal processing. Section |V] 
concludes this work. 
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FIG. 1. Examples of on-site (a-c) and off-site (d-f) gap-soliton 
solutions and the corresponding potential V in the case of the 
local self-focusing Kerr nonlinearity, depicted by black and 
gray lines, respectively. The propagation constant of each 
mode is indicated by labels A through F in panel (g), where 
the relations between the propagation constant and power 
are plotted by the solid and dashed curves, respectively, for 
the on-site and off-site modes. Panel (i) shows the blow-up 
of the area enclosed by the dashed box in (g). Here and in 
other figures, shaded areas represent Bloch bands bordering 
the first finite bandgap. 



II. THE GAP-SOLITON FAMILY WITH THE 
LOCAL SELF-FOCUSING NONLINEARITY 

Considering a wave packet propagating along the rj axis 
in the nonlinear PhC structure, we assume that the em- 
bedded medium gives rise to two kinds of the nonlinearity 
simultaneously. This system is modeled by the modified 
nonlinear Schrodinger equation jisj . 

= + ^(0* + + p\^\H, (f ) 

--rf|? = l*P- (2) 

where ^'(?/, £,) is the slowly varying amplitude of the elec- 
tric field, V{^) is the periodic potential, and a pertur- 
bation of the refractive index, n(^), accounts for the 
diffusive nonlinear response with the nonlocality degree 
(which scales as a squared nonlocality radius) designated 
by parameter d. Sign parameters a, p = +1 and —1 cor- 
respond to the self-defocusing and self-focusing nonlin- 
earities, respectively. Below, we fix cr = -t-1 and p = —1 
for a system with the nonlocal self-defocusing and local 
self-focusing nonlinearitics. 

Stationary solutions with propagation constant —p, are 
looked for as 5'(?7,C) — $('^)e~''"'. For periodic po- 
tential V{S,) = 4sin^(fcoa;) with ko = 1 fixed by rescal- 
ing, the linearized version of Eq. ([T]) gives rise to the 
bandgap structure, with the first finite bandgap being 
f .3047 < p < 3.f896. In this bandgap, GSs can exist in 
form of on-site and ofF-site modes being supported solely 



by the local nonlinearity, i.e., with ct = |i9|]. For the 
self-focusing Kerr nonlinearity, we display examples of 
both the on-site and off-site GS solutions in Fig. [T] The 
relationships between p and the power of these modes, 
P — c?^|$p, are shown in the form of the bifurcation 
curves in Fig. [Ug), with labels A through F referring to 
typical GS mode profiles in subplots (a-f). 

Under the self-focusing nonlinearity, these GSs bi- 
furcate from the upper edge of the first bandgap (in 
other words, from the lower edge of the second finite 
Bloch band), where the spatial-dispersion law features 
the anomalous sign [41i] . In the entire first bandgap, 
the odd on-site GS modes, shown in Fig. [IJ a-c), feature 
two major peaks within one lattice cell . Tails of these 
modes become conspicuous when the propagation con- 
stant moves close to the upper edge of the bandgap. On 
the other hand, a shape transition is demonstrated by 
the off-site even GS modes. Near the lower edge of the 
bandgap, the GS solution has a single major peak coin- 
ciding with a local maximum of the periodic potential, as 
shown in Fig.[T](d). By tracing the variation of the prop- 
agation constant and power along the relation for these 
off-site modes shown by the dashed line of Fig. [IJg) to 
point F, it is seen in Fig.[T](f) that the center of the modal 
profile breaks into two peaks. The double-peaked solu- 
tion is formed due to the balance between the repulsive 
potential barrier and self-trapping induced by the Kerr 
nonlinearity, as the characteristic self-trapping length be- 
comes larger than a half of the lattice period (roughly 
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equivalent to the width of the barrier potential). In be- 
tween, we observe a smooth shape transition of the GS 
from the single-humped shape in Fig. [ijd) to a nearly 
fiat-top one in Fig.[T][e), and, finally to the double-peaked 
mode in Fig. [IJf). 

The off-site GSs exist when the nonlinear self-trapping 
is stronger than the repulsion induced by the potential 
barrier, giving rise to an effective a potential well holding 
the localized modes, as seen in Figs.[ljd-f). On the other 
hand, the effect of the lattice potential is stronger than 
that of the nonlinearity in the case of the on-site modes, 
which, together with the contribution of the gradient en- 
ergy, determines their shapes in Fig.[lja-c). When the re- 
fractive index correction gets weaker in accordance with 
the power reduction, the concave lattice potential create 
a potential barrier to tailor and split the wave function 
into a form of the states that similar to a binding profile 
from two on-site modes in Fig. [IJb). 



III. ON-SITE GAP SOLITONS UNDER 
COMPETING LOCAL AND NONLOCAL 
NONLINEARITIES 




FIG. 2. Examples of on-site (a-e) gap solitons in the case of 
the competing local and nonlocal nonlinearities (black lines) 
and the corresponding potential V (gray lines), for different 
degrees of the non-locality d — 0.1, 1, 40, and the case of p = 
(local only). The corresponding values of propagation con- 
stant /i are labeled (A-E) in panel (f). The dashed line in (f) 
is the asymptotic curve for d = 0.1. 



In this section, we introduce the self-defocusing non- 
local nonlinearity, setting u = 1 in Eq. ([1]). Obviously, 
in this case the total nonlinearity cancels out to zero in 
the limit of d = 0. At d > 0, the overall nonlinearity is 
a self-focusing, because the diffusive nonlocal kernel pro- 
duces a spatially wider and less intensive perturbation of 
the nonlinear refractive index, in comparison with that 
corresponding to the local nonlinearity. Then, similar to 
the situation in the linear model outlined above, one may 
expect the corresponding GS modes to bifurcate from the 
upper edge of the first bandgap, where the effective spa- 
tial dispersion is anomalous. In Fig. [21 we show examples 
of on-site GSs supported by the competing local focusing 
and nonlocal defocusing nonlinearities for different de- 
grees of non-locality, d = 0.1, 1, 40, and the case of p = 
(the local nonlinearity only)for the comparison. 

For different degrees of the nonlocality, the profiles of 
the on-site GS solutions vary slightly, remaining similar 
to their counterparts in the local model. On the other 
hand. Fig. [2Jf ) demonstrates that the power required for 
the formation of the GSs is higher in the model with 
the competing nonlinearities than in the local one be- 
cause the net nonlinearity is effectively reduced by the 
competition of the nonlocal nonlinear response with the 
local Kerr term. In the limit of the ultimate nonlocal- 
ity, d — >■ oo, the -P(/i) curve converges to that in the local 
model (with p = 0) . The latter feature is explained by the 
fact that, in this limit (which corresponds to the model of 
the so-called "accessible solitons" [16]), the nonlocal non- 
linear response amounts to a weak constant background, 
which shifts the propagation constant by the amount pro- 
portional to P/Vd [15] . 

In Fig. [H[f), the bifurcation curve shows an abrupt 
change of the slope both for d = 1 and d = 0.1 (in the 




FIG. 3. Amplitude growth rate of small-perturbation modes 
upon the gap solitons in the case of the competing local and 
nonlocal nonlinearities for the degrees of the non-locality d — 
0.1. The inset shows the blow-up of stability analysis along 
with the corresponding power-propagation constant relations. 

latter case, near point E). In a small region to the left of 
this point, the slope changes its sign to dfi/dP > 0, and 
the corresponding on-site GSs are unstable. It is known 
that the bifurcation curves for GS families supported 
by the local self-focusing nonlinearity obeys an inverted 
" anti- Vakhitov-Kolokolov" criterion Q, dfi/dP < 0, un- 
der which the on-site GSs are stable (see also Ref. [30{). 
To elucidate, we analyzed the stability of the numeri- 
cally found GS families by considering small-perturbation 
modes and calculating their eigenvalues. In specific, we 
show in Fig. [3] the linear stability analysis spectrum over 
the entire first band gap regime. It is clearly illustrated 
in the inset of Fig. [31 unstable small-perturbation eigen- 
modes are obtained only within the region where the 
power dependence d/x/dP > and therefore confirms the 
inverted " anti- Vakhitov-Kolokolov" criterion. The evo- 
lution of these on-site gap solitons generated by direct 
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numerical simulations further shows the collapse of on- 
site solitons falling out of the " anti-Vakhitov-Kolokolov" 
regime. 

When the propagation constant of the GS is close to 
the upper edge of the bandgap (on the right-hand side of 
the slope-change point), the slope of the becomes 
negative again. In that case, the GS is broad, spanning a 
few lattice periods, and resembles gap wave modes [3lj, 
as seen in Fig.[2fe). The abrupt slope change smooths out 
with the increase of the nonlocality degree, disappearing 
at d ~ 5, due to the fact that the nonlocal perturba- 
tion of the refractive index becomes small for the strong 
nonlocality. 

For a sufficiently small degree of the nonlocality, d. 
Fig. [2] shows that the power required to form the GSs in- 
creases with the decrease of d, diverging at c? — )• 0. In the 
regime of the weak nonlocality, the effective competing 
nonlinearity can be approximated, to the first order in d, 
as [H - l^-p^- K d (I^P)^^, which implies that 

the soliton's power scales as 1/d. An asymptotic curve 
based on this approximation is shown by the dashed line 
in Fig. Eff), to illustrate the bifurcation of the on-site 
GSs near the upper edge of the first bandgap. 
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FIG. 4. (Color online) Off-site gap solitons in the model with 
the competing local self-focusing and nonlocal self-defocusing 
nonlinearities for the points marked in panel (c) : A and B (a) , 
C and D (b). The bifurcation curves for two very different 
values of the degree of the nonlocality, d = 0.05 and 20, are 
shown in (c) to indicate the cases with and without the new 
branch. Thin and thick portions of -P(m) curves pertain to the 
solitons with single- and double-peak profiles, respectively. 



IV. OFF-SITE GAP SOLITONS UNDER 
COMPETING LOCAL AND NONLOCAL 
NONLINEARITIES 

Next, we aim to study off-site GSs in the first finite 
bandgap under the action of the competing nonlineari- 
ties. As mentioned in Section [TTl the profile of the off- 
site GS solutions changes from single-humped to double- 
humped as propagation constant /i approaches the upper 
edge of the bandgap. To indicate the change of the pro- 
file caused by the introduction of the competing nonlocal 
nonlinearity, in Fig. |4jc) we use thin and thick lines to 
distinguish portions of the P(/i) curves representing such 
single- and double-humped profiles. Again, in the limit 
of the strong nonlocality, the correction to the refrac- 
tive index induced by the nonlocal nonlinearity is widely 
spread in the space and very small with respect to the 
effect of the Kerr nonlinearity, which makes the compe- 
tition negligible. For example, for d = 20 [see Fig. HJc)], 
both the P(/i) curve and the corresponding shape transi- 
tion of the off-site GSs are close to their counterparts in 
the local model. A smooth transition of the GS solutions 
from the single-peaked shape to a fiat-top one, and then 
to the double-peaked (in-phase) shape can be traced, in 
the latter case. 

However, for a smaller degree of the nonlocality, which 
makes the self-defocusing nonlocal response comparable 
to the Kerr nonlinearity, the analysis reveals the exis- 
tence of more than one branches of the off-site GSs. For 
small values of d, such as d = 0.05 shown in Fig.|4l[c), one 
branch (the bold red curve in the figure) extends continu- 
ously from the upper to lower edge of the bandgap. The 
modal profile for this branch remains double-humped. 



such as the one shown by dashed lines of Fig. HJa) and 
(b) for markers B and D in (c), respectively. Besides this 
branch of the double-humped GS modes, there is a sepa- 
rate branch representing solutions with a lower formation 
power and single-humped profile, as shown by solid lines 
in Fig. |D(a) and (b) for markers A and C in panel (c) . 

Thus, the single- and double-peaked GSs, which consti- 
tute parts of the single GS family in the local model, split 
into two disjoint families in the model with the appre- 
ciable competition between the local and nonlocal non- 
linearities. In the intermediate case, corresponding to 
a moderate degree of the nonlocality, such as for d = 2 
shown in Fig.lS^b), the P{n) curve for the double-humped 
GSs break into two segments: one starts from the lower 
edge of the bandgap and ends at a cutoff point corre- 
sponding to marker B in the figure, and the other starts 
from the upper edge and ends at another cutoff point, 
which is designated by marker A. The cutoff points 
may be accounted for by bifurcations involving additional 
higher-order modes, which we did not aim to find in this 
work dealing with fundamental on- and off-site GSs. An- 
other P(/i) branch chiefly represents the single- humped 
GSs, but it also contains a portion to the right of point 
A, which corresponds to double-humped modes. 

The inset of Fig.O^b) clearly shows that there are three 
branches of double-humped modes. The first one bifur- 
cates from the upper edge of the bandgap and ends at 
the point marked A; the second one extends from the 
lower edge of the bandgap and terminates at point la- 
beled by V; the final branch bifurcates as the flat-top 
solution from the point marked the asterisk (*) and ter- 
minates at the site marked □. The first two branches 
that bifurcate from either edge of the bandgap abruptly 
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terminate inside the bandgap, where the characteristic 
width of the nonhnear response is larger than or compa- 
rable to half the lattice period, which makes the balance 
between the nonlinear and lattice-induced effects impos- 
sible. The last branch, which starts as the flat-top mode, 
ends due to the divergence of the total power as a result 
of vanishing nonlinearity, similar to cutoff considered in 
Ref. [42]. 

Even though off-site GSs are in general unstable both 
in local ^9^ and nonlocal [Ts'l nonlinearities despite the 
inverted " anti-Vakhitov-Kolokolov" criterion, the corre- 
sponding instability growth rate is proportional to the 
GS's power after a certain threshold value [13, [il. Due 
to the unstable nature of off-site GSs, we study the in- 
stability of GSs by the linear stability analysis for the 
branches off-site GSs and identify the final state of these 
off-site GSs by beam propagation simulation. The linear 
stability spectrum in Fig. [5] shows the amplitude growth 
rate of the small-perturbation eigen-modes found upon 
those off-site GS's revealed previously. For a smaller de- 
gree of non-locality, d = 0.05, illustrated in Fig. IUJa) the 
two distinct branches are unstable and the corresponding 
eigen-modes upon single-peaked GSs have higher ampli- 
tude growth rate which diverges near its cutoff point in 
the band gap region. The GSs on these two branches 
collapse fast as they propagates due to a larger growth 
rate. Likewise, in Fig. injb) when a moderate degree of 
non- locality, d = 2, is considered, the single-peaked GSs 
branch acquires a stronger instability than that of the two 
separated double-peaked GSs branches. The two double- 
peaked GSs branches, though both feature the worst in- 
stability near their cutoff points in the band gap region, 
reflect very different relations to the corresponding GS's 
power, which we believe is stemmed from the cutoff of 
the branches. 

The beam propagation simulations further illustrate 
that GSs of the single-peaked branch and double-peaked 
modes with a higher energy fall into collapsed states. 
Nevertheless, mode conversions from unstable off-site 
GSs into stable on-site GSs are observed for the branch 
plotted in blue in Fig. [SJb). Three examples of beam 
propagation simulations are illustrated, resulting in ei- 
ther collapsed states Fig. [7] (a) and (b) or a mode con- 
version behavior Fig. [7] (c). Information such as mode 
transition or conversion is beyond what linear stability 
analysis may reveal. Even though it is also believed that 
interesting dynamical behavior associated with the GSs 
can be delineated through a direct beam propagation 
simulation, yet to be more focused, a thorough inves- 
tigation of propagation behavior goes beyond the scope 
of this work. 

The nonlocal nonlinearity competing with the local 
Kerr term not only reduces the strength of the nonlin- 
earity but also tailors the effective nonlinear response to 
induce binding forces outside of the lattice-potential bar- 
riers, which is a more favorable environment for the exis- 
tence of double-peaked modes. Therefore, the three GS 
branches outlined above are formed owing to the inter- 
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FIG. 5. (Color online) (a) Off-site GSs under the competing 
local self-focusing and nonlocal self-defocusing nonlinearities, 
corresponding to the points A ,B and C in panel (b), where 
the P(^) curves for the nonlocality degree d = 2 are shown. 
The thin and thick lines distinguish portions corresponding 
to single- and double-peaked modes, respectively. Inset in (b) 
is a blow-up of the region around the marked points. 
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FIG. 6. (Color online) Eigenvalues of small-perturbation 
mode upon the Off-site GSs under the competing local self- 
focusing and nonlocal self-defocusing nonlinearities, of which 
the degree of non-localities are (a) d — 0.05 and (b) d — 2. 
The color and line style to distinguish each branch is as is 
defined in the corresponding P(/j.) curves in in Fig. [H^a) and 
Fig. [5jb), respectively. 



play of the potential barrier created by the lattice poten- 
tial and the binding potential induced by the effectively 
reduced Kerr nonlinearity. In this case, the GS solutions 
belonging to the branch originating from the lower edge 
of the bandgap have a larger amplitude and are tighter 
localized. The width of the corresponding response range 
of the self-defocusing nonlocal nonlinearity is larger than 
that of the Kerr response. Such a double-peaked branch 
cannot exist in the strong-nonlocality limit, and we nu- 
merically find that values of the nonlinearity degree sup- 
porting this branch are bounded by d < 9.2. As the 
power decreases (the propagation constant increases) , the 
widening of the the GS mode makes the overall nonlin- 
ear response effectively local, suppressing the capability 
of the nonlocal nonlinearity to tailor its response to a 
shape necessary for supporting the solitons. Then, when 
the GS width becomes smaller than or comparable to half 
the lattice period, the double-peaked modes cease to ex- 
ist because the balance between the lattice potential and 
nonlinearity-induced perturbation of the refractive index 
supports only single-peaked modes. 

To present a clear description of shape transitions for 
the off-site GSs, we replot the relationship of power P 
versus the nonlocality degree, d, for a fixed propagation 
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FIG. 7. (Color online) Three examples of the beam propaga- 
tion simulations for off-site GS solutions shown in Fig. [S][a), 
in which modes A, B and C correspond to intensity plots in 
(a), (b), and (c), respectively. 



constant /i in Fig.|Sl For a smaller propagation constant, 
such as /i = 1.5 in Fig. [5Ja), the power necessary to sup- 
port a single-humped ofF-site GS is always lower than 
that of its double-humped counterpart, for all values of 
d. Moreover, above a critical degree of the nonlocality, 
d = 7.3 in this case, only a single-humped GS can be 
found. Moving fi into the center of the bandgap - for in- 
stance, taking fi — 2.69 in Fig. Hljb) - the critical degree 
of the nonlocality reduces to c? = 1.9, and, above an- 
other critical value, d = 2.35, the single-humped mode 
transforms into a double-humped one, as in the local 
model. Increasing the value of the propagation constant 
to fj, = 2.721855, the two P{d) curves merge at d = 0.3 
in in Fig. [Sfc). For a larger value of the propagation 
constant, such as /i = 2.75 in Fig. [H^d), the two curves 
intersect at a critical value d — 0.09763. Above this crit- 
ical point, the power for the double- humped off-site GS 
becomes lower than for a single-humped one. In this case, 
the profile of the off-site GS mode can be switched from 
single-humped into double-humped by adjusting the non- 
locality degree, d. 

To facilitate the understanding of the present picture, 
we consider the plane of the propagation constant fi and 
power P for the GS solutions in Fig. [3J varying the 
nonlocality degree d. We start by tracing the evolu- 
tion of the point of the transition from single-peaked to 
the double-peaked shape, marker by the asterisk (*) in 
Fig. ^c) . In the absence of the competing nonlocal non- 
linearity, i.e., at p = in Eq. ([T]), the transition point is 
(/i = 2.145, P = 5.0460), labeled by A in Fig. El which 
also corresponds to the limit of d — !■ cx3. As the non- 
locality degree drops to a critical value, c? = 1 at point 
C {fi — 2.7133, P = 8.8184), the transition point ceases 
to exist (i.e., only sharply peaked single- humped modes 
are supported by the system), merging into to the end 
point of the doubled-humped-mode branch, marked by 
□ in Fig.O Increasing the nonlocality degree from d = 1 
at point C, the end point □ in Fig. [5] merges into the 
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FIG. 8. (Color online) The power, P, versus the degree of the 
nonlocality, d, at fixed values of the propagation constant: 
H = 1.5 (a), 2.69 (b), 2.721855 (c), and 2.75 (d), respectively. 
The inset in (c) is the blow-up of the region in the vicinity of 
the merger point. 



other end point A in Fig. [5] at d = 2.127. The latter 
merger happens at point B in Fig. [9l This is the end 
point of the family of the double-peaked modes [an ex- 
ample corresponds to the point marked by A in the inset 
of Fig.[5ljb)] which originates from point B in Fig.|9]at the 
critical nonlocality d — 2.127, and extends toward point 
D, which corresponds to d = 1.86911, where it merges 
into a new branch of the double-peaked modes emerging 
(as long as c? < 9.2) from the lower edge of the bandgap 
at point {fi = 1.3047, P = 19.3262). 



V. CONCLUSION 

In this work, we aimed to study GS (gap-soliton) solu- 
tions in the first finite bandgap of the periodic potential, 
with the nonlinearity represented by the competing lo- 
cal self-focusing and nonlocal self-defocusing terms. The 
two terms are balanced so that, in the limit of the zero 
nonlocality radius, they exactly cancel each other. While 
keeping the effective interaction self-attractive, the exis- 
tence, stability, and bifurcation for on-site and off-site 
modes were analyzed numerically. Due to the opposite 
signs of the local and nonlocal nonlinearities, an increased 
power was required for the formation of both the on- 
site and off-site GSs. The competing nonlinearities in- 
duce a region where stable on-site modes obeying the 
" anti-Vakhitov-Kolokolov criterion" near the upper edge 
of the bandgap region become unstable. For unstable off- 
site GS modes, which remain unstable under competing 
local and nonlocal nonlinearities, a complex bifurcation 
pattern with cutoff points was found and explained in 
terms of the transitions between the single-humped, flat- 
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top, and double-humped shapes. By tracing the evolu- 
tion with the change of the nonlocality degree, we have 
shown that it is possible to switch different off-site GS 
modes by manipulating the nonlocal interaction against 
the local Kerr nonlinearity. 



FIG. 9. (Color online) The diagram in the plane of the prop- 
agation constant (/i) and power (P), for the gap sohtons, as 
different values of the nonlocality degree, d. The inset is a 
blowup of the gray region. Route AC corresponds to the 
transition point marked by * in Fig. 14] routes BC and BDE 
correspond to the double-peaked modes marked by A and V 
in Fig. [S] respectively; the route above point C corresponds 
to the single-peaked solution. The marked points are A with 
p = 0, and B, C, D and E with d = 2.127, 1, 1.86911, and 9.2, 
severally. 
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